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We present results of a study of _B_r- using tadpole improved gauge-invariant staggered operators. Using three 
ensembles of 16 3 x 32 configurations with varying numbers of dynamical flavors, we observe a small dependence 
on Nf. Using 7 quenched ensembles at different values of /3, we extrapolate to a — 0. 



1. Introduction 

The Bk parameter serves to parameterize the 
weak hadronic matrix element responsible for 
K° — K° mixing. Since this mixing gives us the 
only CP violation observed to date, Bk is a cru- 
cial link between the measured quantity e and 
the parameters of the Standard Model. Lattice 
calculations are well suited for the study of B K 
parameter, and it has by now received much at- 
tention. After an early round of calculations 0- 
||, the statistics have now been raised to a level 
which allows one to examine some of the fine 
points of the calculation, such as checks on the 
reliability of one-loop lattice perturbation the- 
ory |Q, the chiral behavior and nondegenerate 
quark masses f|7|, the dependence of Bk on the 
J and the number of dyna- 
In this note we offer more infor- 



lattice spacing [3 
mical flavors 
mation on these latter two points. 



2. Calculational Setup 



as given in table |T[ The quenched configura- 
tions were generated on the Ohio Supercomputer 
Center T3D, while the dynamical configurations 
with two and four flavors of m q a = 0.01 stag- 
gered fermions were generated on the 256-node 
Columbia machine. The parameters were chosen 
so as to make the scales of the lattices exceedingly 
close (and equal to approximately (2 GeV) -1 ), as 
determined from the p-meson mass in chiral limit 
(see Fig. [l] and Ref . || ) . We employ 9 values of 
(degenerate) valence d and s quark masses from 
05. 
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Table 1 

Ensembles for Nf Study 



N f 


/3 


config 


m p a 





6.05 


306 


0.384(5) 


2 


5.7 


83 


0.384(4) 


4 


5.4 


69 


0.391(7) 



m p a vs. quark 



Figure 1. Data and linear fit for 
mass, for three sets of configurations with Nf 
and 4. (See the talk by D. Chen [9].) 



=0,2 



For the continuum limit study we generated 7 
For the dynamical fcrmion comparison we use ensembles of quenched configurations as listed in 

lattices of geometry 16 3 x 32, with parameters table [| and used 7 to 9 values of m q . 
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Table 2 

Quenched Ensembles for Continuum Extrapola- 
tion 



a 


Geometry 


1 Y conng 


m q 




5.70 


16 3 x 32 


259 


.01 to 


08 


5.80 


16 3 x 32 


200 


.01 to 


04 


5.90 


16 3 x 32 


200 


.01 to 


04 


6.00 


16 3 x 32 


221 


.01 to 


04 


6.05 


16 3 x 32 


306 


.01 to 


05 


6.10 


24 3 x 32 


60 


.01 to 


04 


6.20 


24 3 x 48 


121 


.005 to 


.035 



For creating kaons (at rest) we use a wall of 
U(l) noise on timeslice t = 0, i.e. complex 
random numbers £g at each space point such 
that (t;z£,g) — Sg^. This is statistically equiva- 
lent to computing a collection of delta-function 
sources. In particular, our wall creates only pseu- 
doscalars. We use a lattice duplicated in the time 
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Figure 2. We use periodic boundary conditions 
in space and time, and the lattice is duplicated in 
time direction. 



direction, with periodic boundary conditions in 
space and time (see Fig. ||). Computing propa- 
gators on the doubled lattice, we obtain forward- 
and backward-going propagators which we use for 
computing Bk- That is, if G^(t) is the it propaga- 
tor on the doubled lattice, then our operator cor- 
relation functions are schematically of the form 
G„{t)G„{t + N t ), where N t = 32 or 48. 



We employ three kinds of operators: Landau 
gauge, gauge invariant, and tadpole improved. 
Landau gauge operators are defined by fixing the 
gauge and omitting explicit links in non-local op- 
erators. For gauge-invariant operators we sup- 
ply the links, averaging over all shortest paths. 
Tadpole-improved operators are gauge-invariant 
operators, but with all links rescaled by u^ A , 
where u = P 1 / 4 , P is the average plaquette, 
and A is the number of links needed to connect 
fermion fields. We opted for tadpole-improved 
operators on all configurations, using the others 
on a subset of configurations for checks. 

The matching between continuum and lattice 
operators is of the form 

or it = fa + ^hij log (^) + c^of , 

where 7^ is the one-loop anomalous dimension 
matrix, and Cjj are finite coefficients, which can 
be sizable. We take these from the calcula- 
tions of Refs. [ ^0|JTT| . For the continuum scheme, 
we choose NDR, quoting results either at scale 
(j, = tt/ol or at fx = 2 GeV. We use the MS 
coupling constant ffjjg, defined as X/ a) = 
P/3g are +0.02461-Q.00704JV/. To check how well 
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Figure 3. Bk with (lower points) and without 
(upper points) one-loop perturbative matching. 
The points are artificially displaced horizontally 
for clarity. 
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the perturbation theory works, we computed all 
three operators on a subset of the Nf = 2 ensem- 
ble, finding that after one-loop corrections are put 
in, the matrix elements agree within our statisti- 
cal error. For the bulk of the calculation we used 
tadpole-improved operators exclusively. 
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Figure 4. Data and fit for Bk vs. m K on the 



quenched ensemble, 
physical kaon mass. 
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Figure 5. Data and fit for Bk vs. 
dynamical ensembles. The dashed line shows the 
fit for the quenched ensemble. 



lie below Nf — 0, consistent with the observation 
by other groups that quenching seems to increase 
Bk slightly (see, e.g. ref. Ji2fl). However, the 
Nf = 4 data turn this picture upside down. Fig. || 
shows our final values for Bk, obtained at the 
physical kaon mass and by extrapolation to the 
chiral limit. We see that the interpolated Nf = 3 
result is a few percent higher than quenched. 



3. Results for Nf Dependence 

Figs. [| and || show the results for Bk on three 
ensembles of configurations. Values at 9 quark 
mass points are fitted to the form expected from 
chiral perturbation theory, Bk — a + bm K + 
cm K In m K . The Nf — 4 and Nf — 2 curves are 
similar in shape, while the quenched curve crosses 
between the other two. While this is perfectly 
allowed, we should also inject a small note of 
caution — our ensembles have the same p-masses, 
but these masses are presumably affected to some 
degree by the finite volume. If this effect is siz- 
able and depends significantly on Nf, our curves 
could shift a little. 

Taking the results at face value, we note that 
the Nf = 2 and Nf = results lie nearly on top of 
each other at the kaon mass, consistent with our 
earlier results |^|. Also, most of the Nf — 2 data 



4. Continuum Extrapolation 

Performing the same analysis on the quenched 
ensembles, we obtain the result shown in figure 
0, where we plot Bk(NDR, fj, = 2 GeV) versus 
the scale as determined from m p . The data are 
well fit by the quadratic form Bk{o) — Bk(cl = 
0) + (aA2) 2 + (aA4) 4 , where the scale of the power 
correction parameters turns out to be typical of 
QCD: A 2 « 660 MeV, A 4 w 650 MeV. Alterna- 
tively, we note that we can avoid making refer- 
ence to the possibly problematic m p by using the 
scaling form 

«(/ 3 )= a o( T ^j) ex P(7ip-) 

where we take g here to be the MS coupling. This 
amounts to shuffling around the a 4 corrections, 
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Figure 6. Final results for Bk at physical kaon 
mass and in the chiral limit, vs. Nf. 



Figure 7. Linear (heavy line) and quadratic fits 
to Bk{cl). 



and in practice tends to straighten the data out. 
That is to say, much of the curvature in figure 
[?] might be ascribed to scaling violations in m p 
itself. To quote a final value we make the conser- 
vative choice of a linear fit to the four points with 
j3 > 6.0, and obtain 



B 



K 



=0,N,=0 



.573 ± .015. 



5. Conclusions 

From the dynamical comparison, we find that 
B K {N f = 3) is (5 ±2)% larger than B K {N f = 0). 
Combining with the a — extrapolation we we 
quote our current central value Bk in the real 
world: 

B K {NDR, fi = 2 GeV, N f = 3, a = 0) = .60 ± .02 

Remaining uncertainties include possible finite- 
size effects in the dynamical ensemble, higher or- 
der perturbative corrections in the matching, and 
higher order chiral (m s — rrid) effects. A study of 
hadronic weak matrix elements relevant for e'/e 
using the same techniques and ensembles is cur- 
rently underway. 
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